Tomography of pairing symmetry from magnetotunneling spectroscopy — a case 

study for quasi-lD organic superconductors 
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We propose that anisotropic p-, d-, or /-wave pairing symmetries can be distinguished from a 
tunneling spectroscopy in the presence of magnetic fields, which is exemplified here for a model 
organic superconductor (TMTSF) 2 X. The shape of the Fermi surface (quasi-one-dimensional in 
this example) affects sensitively the pairing symmetry, which in turn affects the shape (U or V) of 
the gap along with the presence/absence of the zero- bias peak in the tunneling in a subtle manner. 
Yet, an application of a magnetic field enables us to identify the symmetry, which is interpreted as 
an effect of the Doppler shift in Andreev bound states. 



One of the most fascinating features of unconventional 
superconductors with anisotropic pairing symmetries is 
that Andreev reflection at surfaces can occur in a variety 
of ways |2j . Namely, at surfaces and interfaces, inter- 
ference takes place between incident and reflected quasi- 
particles, which, when the pairing is anisotropic, experi- 
ence opposite signs of the pair potential depending on the 
situation, since an anisotropic pairing dictates that the 
BCS gap function has to have node(s). The interference 
then becomes constructive, and we end up with the An- 
dreev bound states (ABS) at the surface, which should 
be observed as a zero-bias conductance peak (ZBCP) in 
tunneling spectroscopy. This is recognized as a clear sig- 
nature of an anisotropic pairing ||, ||, q|. The ZBCP 
is observed in various anisotropic superconductors such 
as the high-rc cuprates || Q, |[ ||, [To), another oxide 
Sr2RuC>4 p~T|, |l2|l, and a heavy fermion system UBei3 [fl3|| . 
In this context, it is intriguing to investigate whether 
the ZBCP can be observed in organic superconductors, 
another class of candidates for anisotropic pairing, such 
as Bcchgaard salts (TMTSF) 2 Y (X = PF 6 , C10 4 , etc.) 

For (TMTSF) 2 X, a spin-triplet pairing has been sug- 
gested from an observation of a large H C 2 Jl7[ and an un- 
changed Knight shift across T c Esl. Anisotropic pairing 
with nodes on the Fermi surface has been suggested from 
an NMR measurement Jl9| , while a thermal conductivity 
measurement has reported the absence of nodes on the 
Fermi surface [ pfj[ . Theoretically, an anisotropic p-wave 
pairing in which the nodes of the pair potential can be 
made to avoid intersecting the (quasi-lD) Fermi surface 
has been proposed in an early stage plj [2^, ^3). The 
triplet pairing itself, however, is puzzling in an electron 
mechanism for superconductivity, since the usual wis- 
dom dictates that a triplet pairing should be an outcome 
of ferromagnetic spin-fluctuation exchange, whereas the 
superconductivity lies right next to a 2fcp SDW in the 
pressure-temperature phase diagram for TMTSF E4] . - 



naively, SDW spin fluctuations should favor spin-singlet 
d- wave- like pairing p5| , |26| , p7| . In order to solve this puz- 
zle, three of the present authors have recently proposed 
p8| that the correct pairing may be a triplet /-wave, 
which can dominate over d- and p-waves if we introduce 
2k~F CDW coexisting with the 2k-p SDW, and a magnetic 
anisotropy in spin fluctuations in this quasi-lD material. 

Motivated from the above mentioned experimental as 
well as theoretical controversies, we ask ourselves: can 
we identify the pairing symmetry (p, d, f) from the tun- 
neling spectroscopy in general, and from the ZBCP in 
particular? Il4j Recently, it has been pointed out that 
these pairing symmetries in (TMTSF^Y can be distin- 
guished from the presence (for p and /) or absence (d) 
of the ABS on the surface |2{J. Subsequently we have 
pointed out |[o) that p and / can still be distinguished 
by looking at the shape of the gap in the surface density 
of states in which the zero-energy peak (ZEP) J|l| resides 

In real materials, however, the question is what will 
happen for a general shape of the Fermi surface. In 
(TMTSF^Y, the warping of the quasi-lD Fermi surface 
is actually controlled by the ratio between the transfer 
within the chains and several types of transfers across 
the chains. So, if the tunneling spectrum is sensitive 
to the warping, identification of the pairing symmetry 
will be marred, so we will have to devise some in situ 
way of probing the symmetry from the tunneling spec- 
troscopy. In the present study, we first show that the 
tunneling spectrum is in fact sensitive to the warping. 
This is due to the degradation of the symmetry in the 
Fermi surface against k a — k a in k space, and the ap- 
pearance/disappearance of the ZEP may even be inverted 
between d and /-waves by a small change in the hopping 
integral. 

We then show that we can overcome this difficulty by 
applying a magnetic field as an in situ control. In a mag- 
netic field screening current affects the ABS spectrum, 
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which is called the Doppler shift, and this can split the 
ZBCP @, We find that ZBCP splits into two 

for d-wave, while it does not for p and /. The way in 
which the Doppler shift occurs reflects the shape of the 
Fermi surface, so the magnetotunncling spectrum pro- 
vides a unique kind of tomography of the gap function 
symmetry on the Fermi surface. 
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FIG. 1: (a) The hopping integrals and the surface normal to 
the o-axis are shown. Cooper pairs between sites separated by 
mi lattice spacings are illustrated for (b) triplet p- or singlet 
d-wave (m p — = 2) and (c) triplet /-wave {rrif — 4). 



Three of the present authors have derived in Ref. |28| 
plausible pairing symmetries by considering pairing in- 
teractions arising from spin and charge fluctuations. For 
discussing the tunneling spectroscopy for a given pair- 
ing symmetry, it is more convenient to start from the 
extended Hubbard model that incorporates effective at- 
tractions for the given pairing symmetry. So we have 

^ = — tij C ia C jc ~ L JjC i<T Cia 

(i,j),<7 ia 

\i-j\=mi,<7,a' 

where cj creates a hole with spin a =|, J, at site i = 
{ia,ib)- V is the effective attraction that is assumed to 
act on a pair of electrons when they are mi lattice spacing 
apart (along a axis when the system is quasi-lD), where 
mi = 2(4) for I = p, d(/)-wave pairing [see Fig. |l|(b)(c)]. 

As for the lattice we take here an anisotropic (4 <C t a ) 
2D system. To warp the shape of the Fermi surface 
we introduce a diagonal transfer t' ', so the system is an 
anisotropic triangular lattice. The sample edge is as- 
sumed to be _L a axis. In what follows we vary t'(< 0) in 
the range < \t' /t a \ < 0.2 with a fixed tb/t a = 0.1, which 
covers three typical shapes of the Fermi surface (Fig. |^). 
In actual (TMTSF) 2 X salts, \t'\ may be similar to or even 



greater than tb for some anion X at low temperatures 
The chemical potential \i is set (/i 
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-1.41t„) 

to make the band quarter-filled as in TMTSF. 

In the latter half of the paper we apply a magnetic 
field parallel to the c-axis (_!_ a, b). Since the penetration 
depth A > 12000 A|7| in (TMTSF),A is much greater 
than the coherence length £ ~ 500 A(= 150a) g, 
the vector potential can be taken to be A = (0, H\, 0) [||. 
Thus the quasiparticlc momentum kb in the b direction 
changes as k' b = kb + H/(tt(;H c ), where H c = (j> / (-k 2 £\) 
(= 3.5 x 10~ 3 T for TMTSF) with <j) Q = h/(2e). 

For the sample with an edge, a mean-field study is 
performed by introducing a site dependent pair potential. 
The symmetry of the pair potential is best represented 
with the d vector as || A?/' - Y,° q ''' V '' • 

with d\- = ^/(^ • Vy)l'a( c 
a z are Pauli matrices and ctq 



T 1 I 

a >), where a x ,a y , 



= 1. We consider only 
d z ij{= FPj) with all the other <!% = = = for 
triplet p-wave state, d°(= Ffj) only for singlet d, and 

d«(= FL) only for triplet / §t). 

By solving the mean-field equation for a region with 
A^l(= 10 3 ) sites in the a direction and one site in the b 
direction, we have obtained the eigenenergy E v . In terms 
of the eigenenergy E v and the wave functions u%, v%, the 
Bogoliubov-de Gennes equation for the surface with the 
magnetic field applied field parallel to the c-axis is given 
by 
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with H. tj (kb) = -J2±ai t a S 3,i±a + 2i fc cos(k' b a)Sij ■ 
The pair potential is determined self-consistently as 



4- = -^5>^*tanh 



k b ,v 
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2k B T 



(3) 



The ABS is probed by the surface density of states 
calculated with the pair potential determined self- 
consistently. In order to compare our theory with scan- 
ning tunneling microscopy (STM) experiments, we as- 
sume that the STM tip is metallic with a constant den- 
sity of states, and that the tunneling occurs only for the 
site nearest to the tip. This has been shown to be valid 
through the study of tunneling conductance of unconven- 
tional superconductors ||. The tunneling conductance 
spectrum is then given at low temperatures by the nor- 
malized surface density of states B, 
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2k B T 

J2 [K\ 2 t(" - E v ) + K\ 2 5(lo + E u )] . (5) 
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Here ps(w) denotes the surface density of states for the 
superconducting state while pn(u)) the bulk density of 
states in the normal state. 
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FIG. 2: The surface density of states at the surface normal 
to the a-axis in H = for (a) d-wave and (b) /-wave. 



Let us first examine the case of zero magnetic field. 
The shape of the gap in the surface density of states dis- 
played in Fig. H is U-shaped for p-wave, while V-shaped 
for d and /. This is because the nodal lines (displayed 
in Fig. U by the dashed lines) in the pair potential in- 
tersect the Fermi surface for d and /, while they do not 
for p-wave, as noted in Ref. ^ for the case of t' = 0. If 
we turn to the ZEP, p and /-waves have peak [^9], |30"| , 
which is due to the fact that incident and reflected quasi- 
particles normal to the surface feel opposite signs of the 
pair potential, which results in a formation of the ABS 
(whereas oblique incidence is required for d-wave). The 
situation does not change when t' is turned on because 
the nodes of the p-wave lies away from the Fermi surface, 
so the warping of the Fermi surface has little effect. 

On the other hand, the situation is not so simple for 
d and /-waves, where the nodes of the pair potential 
intersect the Fermi surface. For t' — 0, the ZEP appears 
for /, and does not for d as mentioned in Refs. and 
|30| . This is because when the Fermi surface is symmetric 
against injected and reflected quasiparticles 

always feel opposite (same) signs of the pair potential 
for /(d)-wave [see Fig. ||(a)]. This time the situation 
does change when we turn on t', which makes the Fermi 
surface asymmetric against k a <-> — k a , so that some of 
the Andreev reflection [the blue area in Fig. ||(b), left 
panel] lead to ZEP also for d-wave as seen in Fig. 0. For 
\t'/ta\ — 0.15 in particular, the warping of the Fermi 
surface is such that injected and reflected quasiparticles 
mostly feel opposite (same) signs for <£(/) [see Fig. 0(c)], 
so the situation for the ZEP is inverted as seen in Fig. ||. 
Consequently, the surface density of states for d- and /- 
wave pairings look similar (having V-shaped gap + ZEP) 
at around \t'/t a \ ~ 0.1, so that it should be difficult to 
distinguish the two by tunneling spectroscopy. 

This is where the in situ control of the tunneling by 
magnetic field comes in. Namely, we now apply magnetic 
fields along the c-axis. In Fig. we show the magnetic 
field dependence of the surface density of states for d and 




FIG. 3: Fermi surface fixed in t b /t a = 0.1: (a) t' = 0, (b) 
t'/ta = -0.08, and (c) t'/t a = -0.15. The d-wave and /- 
wave pair potential are shown in the left and the right panels, 
respectively. + and — denote the sign of the pair potential, 
and the dashed lines represent the nodal line. The blue areas 
indicate the region where the gap function has opposite signs 
at the ends of the dotted lines, so the states at the ends of 
the dotted lines correspond to injected and reflected states 
of the quasiparticles that lead to Andreev bound states and 
thus zero-bias peaks. In (b), the red areas indicate the region 
where VFb has opposite signs for injected and reflected states. 



/-waves for t' = — 0.08t a . The ZEP for the <i-wave is seen 
to split with the magnetic field, while that for / (or p; 
not shown) does not. The magnitude of applied magnetic 
field here is H < 30H C , where 30H C roughly corresponds 
to H^ 2 (^ 0.1 T), i.e., the upper critical field parallel to 
c-axis, of (TMTSF) 2 PF 6 0. 

The result can be interpreted as follows. For d-wave, 
the ZEP is mainly formed by the quasiparticles having 
kb tt/2, where fFb^a) and i^f&( — k a ) have the same 
sign [see Fig. ||(b), left panel]. This situation is simi- 
lar to those studied previously ||, |[ |||, where the 
ZEP is strongly affected by the Doppler shift. Namely, 
the magnetic field gives in this case the injected and 
reflected quasiparticles additional phase shifts with the 
same sign, which degrades the constructive interference 
for the formation of ZEP. By contrast, for p and /-waves 
quasiparticles with kb — mainly contribute to the for- 
mation of the zero-energy states jij |Il) . For kb = 0, 
VFb(k a ) and Vpb(— k a ) have opposite signs regardless of 
the value of t' [see Fig. ||(b), right panel]. In such a 
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case, the magnetic field gives the injected and reflected 
quasiparticles additional phase shifts with opposite signs, 
which almost cancel out when added, so that the for- 
mation of ZEP is barely affected |l2| . This is physically 
why we can distinguish d and /-waves through the ap- 
pearance/disappearance of ZEP splitting in the presence 
of magnetic field. 
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FIG. 4: The surface density of states normal to a axis for 
\t' /t a \ = 0.08 in the presence of the magnetic field for the (a) 
d-wave and (b) /-wave. 



Although we have here exemplified the tunneling spec- 
tra for the quasi-lD model, the physics involved (i.e., the 
signs of the pair potential and the signs of vp for the 
incident and reflected particles) is quite general, so that 
the principle for distinguishing various singlet and triplet 
pairing symmetries should apply for other, non-cubic (or 
non-tetragonal) systems. 
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